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We investigate the decay process from a time dependent 
potential well in the semiclassical regime. The classical dy- 
namics is chaotic and the decay rate shows an irregular behav- 
ior as a function of the system parameters. By studying the 
weak-chaos regime we are able to connect the decay irregu- 
larities to the presence of nonlinear resonances in the classical 
phase space. A quantitative analytical prediction which ac- 
counts for the numerical results is obtained. 
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I. INTRODUCTION 

In the last years, many conjectures has been put for- 
ward, and tested in various system models, to answer 
the fundamental question of the "quantum chaos prob- 
lem" : what is the signature of classical chaos in the quan- 
tum world? Among these, one of the most intriguing is 
the idea according to which the classical chaos can in- 
duce large scale fluctuations on a genuine quantum phe- 
nomenon such as the tunneling process. Starting from 
the seminal paper of Davis and Heller ||I| , who first noted 
the occurrence of coherent tunneling between regular tori 
separated by a chaotic region, the influence of classical 
chaos on quantum tunneling has been verified in many 
systems and is now accepted in the literature as a fin- 
gerprint of classical non integrability. It is very simple 
to describe this effect. Let us consider a system which is 
classically chaotic and invariant under a symmetry oper- 
ation, like for example the space inversion. If the classi- 
cal system supports a regular torus, by symmetry, there 
might be also a second torus which is distinct from its 
symmetric partner, like for instance two symmetric tori 
encircling the bottom of the two wells of a double-well 
potential. Moreover, let us suppose that the two tori are 
large enough to support quantum states. In this condi- 
tion, the quantum system will show coherent tunneling 
between the states located in the two symmetric tori. If 
now one system parameter is changed (e.g. h), contrary 
to the expectations of ordinary semiclassical analysis, the 
tunneling rate shows strong irregularities which can in- 
crease or decrease the rate by orders of magnitude. This 
effect does not show up in a system whose classical coun- 
terpart is integrable. 

The tunneling fluctuation is usually interpreted in 
terms of assisted processes, or, using a widespread ter- 
minology, as "Chaos Assisted Tunnehng" (CAT) [|-0. 



An intuitive view of the CAT process is as follows. The 
presence or regular and stochastic motion in the classi- 
cal phase space corresponds, from a quantum point of 
view, to the possibility of having two kind of states: reg- 
ular ones localized inside the symmetric tori and chaotic 
states, which being extended through the chaotic region 
have a not negligible overlap with regular regions. The 
fluctuations in the tunneling rate are thus explained in 
terms of a three state tunneling process. The quantum 
particle first tunnels from the localized state to an ex- 
tended chaotic one and then from this to the state located 
in the symmetric torus. 
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FIG. 1. Sketch of the typical behavior of the energy lev- 
els of a classically chaotic quantum system as a parameter 
p is changed. In a) the two solid lines describe a couple of 
quasi-degenerate levels of different symmetry. The dashed 
line describes a colliding third level. In b) we show the split- 
ting of the two regular levels. All the units are arbitrary. 

This interpretation is confirmed by the level dynamics 
of the tunneling system. A typical situation is sketched 
in Fig. |l| where we can observe the change of two quasi- 
degenerate levels, which correspond to the pair of tun- 
neling regular states, as a system parameter is varied. In 



almost the entire parameter range the sphtting between 
the two states, and so the direct tunneUng probabihty, 
changes smoothly. However, it may occur that, once the 
parameter is changed, a third level (dashed line) crosses 
the two quasi-degenerate levels. In the non-integrable 
case states belonging to the same symmetry class do not 
cross each other. Therefore, the appearance of a third 
colliding states gives rise to an avoided crossing with the 
state of the doublet which belongs to the same symmetry 
class. 

The avoided crossing has a twofold consequence on the 
tunneling process under study. First, in the nearby of 
the crossing we cannot consider the tunneling as a pro- 
cess involving only the two quasi-degenerate states. In 
this condition the standard two state tunneling becomes 
a resonant three state process. Second, since the col- 
liding third state modifies the energy level of only one of 
the doublet states, the splitting of the two levels changes. 
Therefore, given the small value of the energy splitting in 
the semiclassical regime, the avoided crossing produces a 
dramatic modification of the levels splitting and conse- 
quently of the tunneling rate. It is important to point out 
that the rate can increase by several orders of magnitude 
as well as vanish (see the arrow in Fig. n|) according to 
the value of the parameters. Moreover, due to the fact 
that the energy spectrum, in the nonlinear case, does not 
show any regularity, the crossings with a third level do 
not follow a regular pattern, and the overall behavior of 
the tunneling rate appears to be an irregular sequence of 
peaks 0-01 instead of the smooth behavior expected in 
the regular systems. 

These features motivated the widespread idea that 
classical chaotic trajectories can have an active role in 
the quantum process, helping or "assisting" the quan- 
tum particle to tunnel between the symmetric tori. Along 
this line, path integral techniques have been recently used 
to calculate the contribution to the tunneling stemming 
from complex orbits which connect the symmetric regular 
tori through the classical stochastic layer . 

However, a real quantitative theory of the CAT process 
is still lacking. The main reason for this can probably be 
found in the chaotic nature of the third state which pre- 
vents any simple analytical treatment. Moreover, there 
are some aspects of the phenomenon which do not seem 
to fit properly with the intuitive interpretation given by 
the CAT picture. For example, the presence of strong 
decreases in the tunneling rate which, together with the 
enhancements, occur as a result of a parameter change, 
contradicts the idea of a tunnel process being "assisted" 
by chaos. Another controversial aspect of the effect is 
the nature of the third state which crosses the tunnel- 
ing doublet. The key point of the CAT picture is that 
the perturbation in the energy splitting is relevant only 
for those crossings involving colliding third states which 
are located in the chaotic region. However, recently, it 
has been shown that similar strong fluctuations can be 
obtained also in almost integrable systems 0. In this 
condition, the third state responsible for the fluctuation 



is by no means chaotic. It is regular and its crossing with 
the tunneling doublet can be directly related to the clas- 
sical phase-space structure and in particular to the de- 
struction of the regular tori which are transformed into 
classical non-linear resonances. In particular, a fluctu- 
ation (avoided crossing) occurs when the energy of the 
tunneling state corresponds classically to the energy of a 
nonlinear resonance. In other words, in this regime, the 
fluctuations (avoided crossings) are the quantum mani- 
festations of the classical "small denominators" problem. 
This discovery also led to a first analytical prediction 
about the positions of the tunneling fluctuations which 
we shall review in Section IV. 



In this paper we want to assess whether this picture 
applies also to a different tunneling process, namely the 
quantum escape of a particle which has been initially lo- 
cated inside a potential well. From a classical point of 
view it is clear that the particle can overcome the poten- 
tial barrier of the well only if its energy is larger than the 
barrier height, while in the quantum framework the tun- 
neling across the classically forbidden region is always 
present. Clearly this situation is going to be modified 
when, including the ingredient of chaos, we perturb the 
system by adding a forcing term, i.e., a time dependent 
external force. The perturbation disturbs the regular mo- 
tion of the classical particle and, by increasing its energy, 
makes it possible to the particle to escape over the bar- 
rier. In the meanwhile, also the quantum process of tun- 
neling changes due to the modification of the potential 
and both the processes contribute to the decay 1 11|] . Our 
purpose is to choose a region of the system parameters 
where the classical and the quantum contribution to the 
decay can be separated, in order to study the proper- 
ties of the latter process in connection with the chaotic 
features of the classical phase space. 

By studying a decay process, we shall address a infinite 
system, with a continuous spectrum, and this will prevent 
us from using standard methods, like the diagonalization 
of the Floquet dynamic operator [Q , to obtain directly 
the level splittings responsible for the tunneling. Any- 
way we shall be able to analyzes the system by resorting 
to a somewhat simpler method: we shall calculate nu- 
merically the time evolution of a quantum state initially 
located in the potential well and, by studying the decay of 
the population in the well, we shall be able to obtain the 
relative strength of the tunneling as a function of the sys- 
tem parameters. This will allow us to point out the dif- 
ferences between this process of chaos assisted decay and 
the tunneling processes in the presence of chaos. Finally, 
we shall be able to show that the picture which singles out 
the classical nonlinear resonances as the main responsi- 
ble for the fluctuations of the tunneling rate, applies also 
in this context, and that the prediction of Ref. H, illus- 
trated in Section IV, remains valid. 



II. THE MODEL: CLASSICAL DYNAMICS 

In order to analyze the influence of classical chaos on 
the quantum process of escape from a potential well, 
we introduce a simple one-dimensional forced system de- 
scribed by the following Hamiltonian: 



V{q,t) ^Vo{l~cos{2q)) + 

e{l - cos{2q - fit)) q 
V{q, t) — otherwise . 



[-7r,7r] 



(1) 

(2) 
(3) 



The particle is located initially inside a potential well 
which has the form of a sinusoidal function extended over 
two periods, as shown in Fig. |^, and then it is forced by 
a time periodic perturbation which is considered to be 
small compared to the static potential well, i.e., e <C Vq. 
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FIG. 2. The unperturbed potential of Eq. |2[ We also show 
the wave function of the two states used as initial condition 
in the numerical calculations. 

We chose a perturbation term which turns system dl^) 
into a double resonance-like Hamiltonian. This choice is 
dictated by the sake of simplicity. It is indeed clear that, 
as long as we limit our analysis to small perturbations, 
the particular form of the external forcing does not af- 
fect the generic features of the decay process we want 
to study. On the other hand, the adoption of Hamilto- 
nian (M presents many benefits. All the relevant infor- 
mation concerning the dynamic properties of our model 
can be derived from the dynamics of a well known system, 
the double resonance Hamiltonian which corresponds to 
eq. (P with periodical boundary conditions p3|J^ . 

The presence of a periodic perturbation in^ipbreaks 
the integrability of the classical Hamiltonian. The most 
important features of this condition is the appearance of 
nonlinear resonances in the phase space together with re- 
gions characterized by extended chaotic motion (stochas- 
tic layer). The relevance of the chaotic motion depends 
on the strength e of the perturbation term, so that the 



system can be more or less chaotic. In Fig. g we show a 
stroboscopic mapping of the dynamics, namely the posi- 
tion in phase space at fixed intervals of time which are 
integer multiples of the forcing term period T — 27r/r2, 
for a generic weak-chaos case. Some nonlinear resonances 
and the stochastic layer around the separatrix are clearly 
visible. 

The nonlinear resonances are the visible consequences 
of the small denominators problem. These are related to 
the secular terms which appear in the perturbative solu- 
tion of the equation of motion of non-integrable systems. 
In the weak-chaos condition their position in the phase 
space can be obtained by considering the effect of the 
time-dependent term as a perturbation on the dynamics 
expressed by the constant Hamiltonian 



p 
Ho = — + Vb(l - cos(2g))|_^<q<7r 



(4) 



The KAM theorem |15| states that, as long as the per- 
turbating term can be considered small, the main part 
of the phase space remains practically unperturbed and 
that only the tori which are resonant with the forcing 
term are destroyed and replaced by chains of islands like 
the ones shown in Fig. ^ The resonant condition can be 
written as 



TTl 

ujo{E) = —n 

n 



(5) 



where n and m are integer numbers and ujo{E) is the fre- 
quency of the unperturbed motion inside the well which 
depends on the energy E. In our case it is possible to 
express cooiE) in term of the elliptic function K(k) as 

(6) 



uJoiE) =. Tr^Vo/K{k) 
k={E + K))/2Fo . 

The result of Eq. (ph would actually indicate that all the 
tori are destroyed by the perturbation, being the rational 
numbers dense among the real ones. Nevertheless, the 
KAM theorem assures that the effects of the perturbation 
becomes smaller and smaller with increasing the order of 
the resonance, i.e., with increasing the numerator m. 

0.5 

0.4 

0.3 

0.2 

0.1 
p. 
-0.1 
-0.2 
-0.3 
-0.4 
-0.5 






-1.5 



-1 



-0.5 



0.5 



1 



1.5 



FIG. 3. The classical dynamics inside the well. Strobo- 
scopic Poincare map. The values of the parameters are 
Vo = 0.048, e = 0.005, Q = 2. 



This is clearly visible in Fig. ^ where we can only rec- 
ognize the chains corresponding to m = 1, i.e., the 1/5, 
1/6 and 1/7 resonances. However, even for small per- 
turbation, in the neighborhood of the separatrix of the 
unperturbed system the motion is always dominated by 
chaotic dynamics. In other words, trajectories which in 
the absence of perturbation are bounded inside the po- 
tential well, can now overcome the energy barrier and 
eventually escape from the well. In the generic weak 
chaos condition, the stochastic layer around the separa- 
trix is dynamically separated from the phase space region 
corresponding to bounded trajectories by unbroken tori, 
so that the process of escape driven by chaos is limited in 
phase space. Therefore, in the classical case, the decay 
of the population of the well is possible only if some par- 
ticles are initially put inside the stochastic region. This 
process have been extensively studied in the last years in 
various classical and quantum models, the most known 
of which is probably the Hydrogen atom in the presence 
of a strong radiation field [^6|jl7]. In that system the dy- 
namic process described above eads to the ionization of 
the atom. However, in this paper we want to focus on the 
connection between chaos and processes which would be 
classically impossible, such as quantum tunneling. For 
that reason we shall analyze the decay of the well pop- 
ulation for a quantum particle initially put in the phase 
space region which corresponds, even in the presence of 
chaos, to bounded motion. In connection to this it is 
worthwhile to remember that in quantum mechanics the 
situation is never as simple. Whatever the initial condi- 
tion is, the wave function cannot be sharply located inside 
a finite region but exhibits smooth decreasing tails which 
extend over the stochastic layer. Therefore, to keep the 
classical chaotic diffusion process as small as possible, we 
shall consider a weak-chaos regime with a small stochas- 
tic layer like in Fig. ^, and in addition to this, we shall 
study the decay of quantum states well localized inside 
the potential well. 



III. THE MODEL: QUANTUM DYNAMICS 

We studied the quantum decay from the well of 
Fig. by integrating numerically the time dependent 
Schrodinger equation associated with Hamiltonian p). 
This can be done using a FFT splitting algorithm p^ and 
using absorbing boundary conditions p9| | as described, 
for example, in Ref. p^ . 

In order to single out the effect of the chaotic pertur- 
bation on the process of quantum decay, it is necessary to 
choose as initial condition a state localized inside the well 
which has an unperturbed dynamics as simple as possi- 
ble. The eigenstates of Hamiltonian (nl) with e = are 
not useful in this context being system (|l|) an open sys- 
tem (continuous spectrum with stationary states which 
do not have finite support inside the well). We thus re- 
sorted to use as initial condition mctastable states which 



have a negligible internal dynamics and a long enough 
unperturbed life-time inside the well. These are the res- 
onances of the potential well of Fig. defined in the quan- 
tum theory of scattering. In addition to this, it is worth- 
while to notice that, since in the parameter range we 
explored the unperturbed decay probability is negligible, 
the adoption of the eigenstates of Hamiltonian (0) sup- 
plemented by periodical boundary conditions leads essen- 
tially to the same results. Therefore, as a first approx- 
imation we can consider the initial states as stationary 
state of the unperturbed system. For the sake of simplic- 
ity in the following we shall refer to them as "eigenstates" 
of the unperturbed Hamiltonian. 




FIG. 4. The time evolution of the well population. The dif- 
ferent curves corresponds to 20 different values of Q included 
between fi = 1.5 and Jl = 2.5. The values of the parameters 
are Vo = 0.048, e = 0.005, h = 0.025 and the initial state is 
the fourth eigenstate of the unperturbed well. The horizontal 
dashed line correspond to the unperturbed e = case. Note 
the fast oscillations in the population decay which reflect the 
oscillation of the forcing term 

In Fig. we show the time evolution of the popula- 
tion inside the well P{t) = J^l/2 |V'(9:^)P'^'/ for different 
values of the forcing frequency Q. This figures has been 
obtained by choosing as initial condition the fourth eigen- 
state, 1 3), of the unperturbed well. It can be easily real- 
ized that the decay probability can be strongly enhanced 
by the forcing term even in the small perturbation regime 
(e — 0.005). Note that in the unperturbed case the pop- 
ulation decay is not visible in the scale of this figure, the 
population remaining practically unchanged in the stud- 
ied interval of time. This shows that the time scale of the 
unperturbed process is much longer than the maximum 
time which we explored numerically. 

To obtain a more quantitative representation of the 
phenomenon wc could define the decay rate as the inverse 
of the time integral of P{t), namely the inverse of the area 
contained under the curves of Fig. 0. This would imply 
a very long numerical simulation, up to a time where the 
population has completely leaked out. However, we are 
not interested in the absolute magnitude of the decay 
rate, but only in its relative strength as a function of 



the system parameters. Thus we can simply calculate 
the time integral of P{t) up to a certain time tmax and 
measure the rate by studying the quantity 



R= 1 



1 



P{t)dt 



(7) 



as a function of the frequency Q of the forcing term. As 
shown in Fig. ^, R shows a sequence of peaks, similar to 
a resonant dependence on the frequency of the perturba- 
tion. We repeated the calculation for two different initial 
conditions, namely by choosing the third and the fourth 
eigenstate of the unperturbed well. The decay rate for 
the third state is smaller, as expected, since this state 
is more deeply located in the potential well, but in both 
cases we found a similar behavior even if the position of 
the peaks and their intensity is not the same. 




FIG. 5. R, calculated with tmax ~ 400, as function of 
the driving frequency Q. The values of the parameters are 
Vo = 0.048, e = 0.005, h = 0.025. The two figures correspond 
to different initial conditions, Fig. a) by choosing as initial 
state the fourth eigenstate of the unperturbed well. Fig. b) by 
choosing the third. The numbers refer to the classical nonlin- 
ear resonances as discussed in Section IV, The insets contain 
enlarged views. 

Except for the presence of the peaks, the figures shows 
that, as a general tendency, the decay increases as the 



forcing frequency decreases. This can be understood us- 
ing arguments based on the classical dynamics of sys- 
tem (0). Indeed, as shown in Fig. 0, the chaotic features 
of the classical phase space (i.e., width of the stochastic 
layer and of the nonlinear resonances) increase by re- 
ducing ri, thus demonstrating that the forcing term be- 
comes more important when fl gets smaller. Eventually, 
for extremely small fl, the stochastic layer becomes so 
wide that, for the chosen initial conditions, the escape 
from the well via the direct classical process becomes the 
dominant process. Since we are rather interested in the 
quantum mechanism of escape from the well, we shall 
not explore the condition of small f2's corresponding to 
strong classical chaos. 

In addition to this, there is a further reason to limit 
the analysis to fi's not too small. Our evaluation of the 
decay rate, as the integral of the surviving population 
in the well over a finite time, is meaningful only if the 
time of integration is much longer than the period of the 
perturbation, and therefore we shall limit ourselves to 
study the range of large Q {2n/n « tmax)- 

However, it is important to point out that the numeri- 
cal results represented in Figs. ^ and || does not account 
for the process of classical escape from the well, even to 
the smallest used value of fi. A direct calculation demon- 
strates indeed that the classical decay is always negligi- 
ble in the considered parameter range of fl. This can be 
easily assessed by numerically integrating system (|^) us- 
ing a classical initial particle distribution mimicking the 
phase-space representation of the initial quantum state 
(see Section ^ for details). Even for the smallest value of 
fl used to obtain the results of Figs. || and ^, the classical 
population in the well does not practically change with 
time. 

Therefore, the peaks and the general tendency of 
Figs. ^ and |^ are genuine quantum effects, but, while 
the latter can be associated to the increased effectiveness 
of the forcing term in (m) with decreasing il, the former 
effect does not have any simple interpretation. The un- 
derstanding of this will be the subject of the next section. 

To conclude the analysis of Fig. ^ we have to assess the 
role of the quantum decay process in the absence of the 
external perturbation (i.e., e = 0). By direct numerical 
integration we observed that in the unperturbed case the 
decay rate is orders of magnitude smaller than that of 
the forced case. Indeed, for e = 0, the population inside 
the well remains equal to one (within the numerical pre- 
cision) even at i = tmax- This confirmed that within the 
observed time, the chosen initial states can indeed be re- 
garded as "stationary states" of the unperturbed system. 
On the other hand, this result, as well as previous obser- 
vations, suggest that the sudden increase of the decay 
rate is connected to the classical chaos in a way which 
is not directly connected to the classical process of bar- 
rier crossing via chaotic diffusion through the stochastic 
layer. In the next section we shall build up a theory to 
explain the presence of peaks. 



IV. A SEMICLASSICAL ANALYSIS 

The results shown if Fig. S resemble the typical CAT 
behavior: when a system parameter is changed the decay 
rate presents an irregular sequence of peaks on a smooth 
background. It is thus natural to look for a connection 
between CAT systems and our model. First of all, we 
notice that, due to the continuous nature of spectrum, 
the level dynamics of our system can not be simply de- 
scribed in terms of avoided crossings. Nevertheless, we 
believe that the phenomenon underlying Fig. |5| retains 
much of the features of the CAT processes. In particu- 
lar we think that the argument introduced in Ref. Q to 
explain the tunneling irregularities of a quasi-integrable 
system, can be effectively used also in this model. In 
the cited reference, a connection between the peaks in 
the tunneling rate and the position of the nonlinear res- 
onances in the classical phase space was found. In this 
section we shall exploit the same idea. 

The analysis of Ref. § follows the line of Refs. [|l]J2| 
and is based on a semiclassical approximation which 
makes use of simple arguments. We shall assume that 
the classical system is only weakly perturbed by the ex- 
ternal perturbation: the size of the chaotic region is con- 
sidered to be small with respect to the portion of the 
classical phase space covered with regular trajectories so 
that the main effect of the perturbation is the appear- 
ance of chains of nonlinear island in the classical phase 
space as shown if Fig. 0. In this condition, the area of 
the regular region, that is the phase-space region which is 
encircled by the last unbroken torus inside of the stochas- 
tic layer, is much larger than h. In this way, the regular 
region can accommodate several quantum states and the 
semiclassical approximation becomes meaningful. 

As we discussed in the Introduction, all the studies 
about the effects of chaos on tunneling or on decay pro- 
cesses were concerned with the role of the stochastic layer 
and the classical transport therein as the main contribu- 
tion to the barrier crossing. This effect is surely present, 
but its contribution is not always the most important, at 
least in a weak- chaos regime. In fact as shown in Ref. Q, 
by studying the tunneling process in a forced double well 
system in a condition on weak chaos, it has been possible 
to clarify the role of until now disregarded actors: the 
nonlinear resonances which are present inside the regular 
region of the well. These pre-chaotic structures cannot in 
fact contribute to the classical transport over the barrier, 
because they are embedded in a regular region of unbro- 
ken tori, but they can perturb the quantum dynamics via 
a indirect process similar to the one discussed in CAT. 
The process of avoided crossing responsible of the tun- 
nel irregularity in CAT can in fact be present also in 
a weak-chaos regime, with the difference that the third 
state which, by crossing, modifies the tunneling, needs 
not to be chaotic. 

This process, which we think to be a special realization 
of the more generic CAT one, in the weak-chaos regime 



can give a contribution to the tunneling rate modifica- 
tion even more important than the one connected to the 
chaotic region of the phase space. Moreover, due to the 
fact that the perturbing third state is not chaotic, it is 
possible to obtain a quantitative prediction of the avoided 
crossings and hence on the positions of the tunneling ir- 
regularities. We shall now review the derivation of this 
prediction which has been only sketched in |9[| . 

Let us note first of all that the Hamiltonian under 
study is time dependent and that this prevents us from 
adopting the Einstcin-Brillouin-Keller (EBK) quantiza- 
tion conditions 123] in their original form. However, fol- 
lowing the work of Breuer and Holthaus [p2|, who, in 
turn, extended the method of the canonical operator 
as developed by Maslov and Fedoriuk [g4|, it is possi- 
ble to adapt the EBK prescriptions to periodically time- 
dependent systems. The generalization of Ref. [g2| is 
based on the prescription of Arnold |23] and leads to 
semiclassical quantization rules for the Floquet quasi- 
energies and quasi-eigenstates Q. This is made possi- 
ble by a suitable extension of the phase space, including 
the time i as a coordinate and adding the corresponding 
conjugate momentum. In the one-dimensional case the 
semiclassical quantization prescriptions read: 



•^ = ^ f pdq = h{n+ -), 



(8) 



71 



En,m ^ ^Tf; (pdq - Hdt) + nnm. 

^ -^72 



The meaning of the symbols adopted in this expres- 
sion can be explained by remembering that in the one- 
dimensional case the extended space is the tridimensional 
phase space {q,p, t}. A regular trajectory is contained on 
a flux tube in this tridimensional space, flux tube which 
repeat itself periodically along the t direction. Thus J is 
the action associated to the quantized trajectory, 71 is a 
close path winding once around the flux tube and lying 
in the plane at a given time t and 72 is a path stretching 
itself out on the surface of the flux tube such that it can 
be continued adopting periodic boundary conditions. In 
other words, the path 72 in the extended tridimensional 
phase space, moves from an initial point lying on the 
plane i = to a flnal point lying on the plane t = T. 
Note also that we can choose to lay the path 71 on the 
Poincare section of the flux tubes. Finally, if we restrict 
ourselves to consider the closed orbits inside the potential 
well, the Maslov index assumes only the value v = 2. 

Worth of a detailed discussion is the structure of the 
quantized energy En^m (for simplicity we shall use the 
terms energy and state as equivalents of quasi- energy and 
quasi-state). While the index n, which values are flxed by 
the first of Eqs. (^ , has the usual role of principal quan- 
tum number, the index m, and the dependence of En^m 
on this one, reflects the periodicity of the time dependent 
term of the Hamiltonian. In fact, an important aspect 
of the Floquet theory is the Brillouin zone structure of 
the energy spectrum: for each physical solution labeled 



by n we have a infinite series of representative labeled by 
the value of m. Naturally all the physical information is 
contained in the first Brillouin zone < En^m < ^^, or 
equivalently, we can say that any solution of Eqs. (|^) can 
be folded back to the first Brillouin zone by an appropri- 
ate choice of m. 

At this point it is important to recall that the earlier 
formalism represents a valid quantization procedure only 
in closed systems, where the notion of energy levels is 
meaningful, and where the phase space is filled with reg- 
ular tori. However, we are investigating the properties of 
states that lie well inside the stability region, and that 
present a small decay probability (see Fig. ^. In this 
condition, we believe that an analysis which, according 
to the prescriptions of CAT, connects the tunneling ir- 
regularities with the presence of avoided level crossings, 
can retain its validity. Let us thus proceed disregarding 
the continuity of the spectrum and looking for the pres- 
ence of level crossings in the unperturbed hamiltonian 
spectrum as a function of the forcing frequency Q. 

Following Ref. |g^ we can look for the level crossings 
by replacing H with Hq , let us in fact recall that we are 
always considering a perturbative approach. The condi- 
tion of crossing between states n and n' in the Brillouin 
zone yields the following equation: 

Ho{h{n + 1/2)) + nnm = Ho{h{n' + 1/2)) + hnm' . 

(9) 

This equation can be simplified if we assume that h is so 
small as to make negligible the quantity h{n — n'). If we 
now expand the RHS of this equation around this small 
parameter, we obtain 



dn 



(n — n) + Mlm — Mlm 



which can be rewritten as 
dHodJ 



dj dn 



{n — n') — hn{m' — m) 



or, by using J = Ti{n + 1/2), 



An 



(10) 



(11^ 



(12) 



^^ is the frequency of the un- 



where loq = a;o(J) — ^j 
perturbed motion as a function of the classical action J, 
An = n' — n and Am = m — m' . The condition of levels 
crossing, in the limit of vanishing h can thus be obtained 
by solving the classical equation which corresponds to 
the condition for the onset of the classical nonlinear res- 
onances. 

This is a very important results, because it means that, 
for sufficiently small h s, there is a correspondence be- 
tween the presence of a nonlinear resonance in the classi- 
cal phase space and a level crossing in the spectrum and 
therefore a correspondence between the tunneling irreg- 
ularities and the nonlinear resonances. In other words, 



the tunneling peaks are the quantum manifestation of 
the non-integrable behavior of the underlying classical 
dynamics. The results of Eq. ( [l^ ) is even more specific, 
in fact is says that the crossings of the two unperturbed 
levels E^ ^ and E^, ^, is related to the superposition be- 
tween the semiclassical quantization torus of one of the 
two states and the nonlinear resonance of appropriate 
order Am/ An. 

This process admits a simple intuitive representation. 
Let for example consider a quantum state located deep 
in the well. Its quantization torus lies inside the well and 
the wavefunction in the semiclassical regime is mainly 
located around this torus. The decay rate in the unper- 
turbed case can be obtained by the usual semiclassical 
calculation of the probability of barrier crossing. If we 
turn the perturbation on in a regime of weak chaos we 
have a high probability that nothing happens, due to 
the fact that most of the phase space inside the well re- 
mains unperturbed. But if we change an external param- 
eter as, for example, the forcing frequency fi, we obtain 
that the nonlinear resonances move in the phase space 
and, for particular values of the parameter, one of them 
can intersect the quantization torus which is destroyed. 
This would be probably reflected in a perturbation of the 
quantum state and thus in a modification of its decay 
rate. This is exactly what is described by Eq. (|l^ ) 

Notice anyway that Eq. (113) involves resonances of any 
order and this implies that the change of a system pa- 
rameter makes a chosen level undergo a virtually infinite 
number of crossings. In other words, when we turn the 
perturbation on, the mathematical curve corresponding 
to a particular energy level becomes extremely complex, 
being fragmented by the avoided-crossing effects in an in- 
finite number of points. However, the perturbation pro- 
duced by the avoided crossings strongly depends on the 
order Am of the resonance and in practice, if we restrict 
ourself to the perturbative regime, the only significant 
crossings are those associated with the first-order reso- 
nances (Am = 1). On the other hand the same happens 
in the classical dynamics, as discussed in Section (y). 

The simple prediction of Eq. ( [l2|) is valid in the strong 
semiclassical limit and is therefore unavailable for a nu- 
merical check, being the smallest value of Ti dictated by 
the computer limitations. The relation between non- 
linear resonances and tunneling peaks has been indeed 
proven in Ref. B by means of a further expansion of 
Eq. (H). To do that we proceed by evaluating the second 
order term in h{n — n'). We obtain: 



On the other hand 



2 dJ An ^ ^ 



dujQ 
Thus we can write 



dajQ dE 
IWdJ 



dujo 



dE 
3|) as follows: 
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(14) 



ujo{E) = 



n 



Ar 



(1 



2 dE ^'V 



o{h^) 



(15) 



Where we can give an analytical expression to ""j^ ' , by 
using Eq. (P), as 



dujQ 
'dE 



1 



4fc2^ 



-(1 



^(fc) 
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(16) 



%K{k)^ k'^K{k)' 

where fc'^ = l — k"^. Eq. ( |l5| ) is the main results of Ref. ||] 
and represents a generalization of the classical nonlinear 
resonance condition, where the frequency ratio is renor- 
malized by means of a quantum correction proportional 
to h. 

By using this prediction we are now able to verify our 
conjecture about the validity of this method also in the 
present case. From Eq. (|l5| ) it is possible to predict the 
position of the peaks of decay. After fixing the values of 
Am Eq. (^5|) can be solved as a function of the energy E 
for several values of n. The graphical solution for the first 
order, Am — 1, crossings is shown in Fig. |6| where the 
thick solid horizontal lines correspond to the semiclas- 
sical energies of the third and fourth eigenstates of the 
unperturbed Hamiltonian (0), i.e., of the states which we 
chose as initial condition in order to obtain the results of 



Section III. The dotted and dashed curves represent the 



quantum renormalized energies of the classical nonlinear 
resonances of different order n (the order is indicated by 
the numbers in the figure). The crossings between the 
horizontal lines and the decreasing curves thus indicate 
the solutions of Eq. (Ilq). Their position should also in- 
dicate the position of the peaks of the decay rate. This 
is in fact approximately true, as one can check by going 
back to Fig. ||, where we indicated the solutions showed 
in Fig. with the numbered vertical dashed lines. 
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FIG. 6. The solutions of Eq. ( |l5| ) represented by the cross- 
ings between the solid horizontal lines, indicating the energy 
of the third and fourth eigenstates, and the energy of various 
nonlinear resonances indicated by their order n (dashed and 
dotted hues). The crossings give the theoretical prediction on 
the position of the peaks of decay and are reported in Fig. M 
as the vertical dashed lines. The values of the parameter are 
Vo = 0.048, e = 0.005, Q = 2, h = 0.025. 



We are now able to justify a posteriori the use of 
the semiclassical theory of this section in the present 
case. The main difference between the physical system of 
Ref. 1^ and Hamiltonian (||) is the energy levels discrete- 
ness. In the present work we cannot speak about level 
crossings and thus the calculations above could look in- 
validated. On the other hand we showed that the avoided 
crossings are nothing more than the trait d'union be- 
tween the nonlinear classical resonances and the peaks of 
the decay rate: the presence of the nonlinear resonance 
produces a level crossing which is reflected in a rate ir- 
regularity. The same happens for the system of Hamil- 
tonian (|^): the classical phase space structures and the 
decay rate modifications are related even if the intermedi- 
ate step is less clear due to the continuity of the quantum 
spectrum. Probably we could find a process similar to the 
avoided crossings, but we do not need to look for it as we 
showed that the quantum-classical connection works. 

The role of the nonlinear resonances in the perturba- 
tion of tunneling seems thus to be established also in 
a system with a continuous spectrum, even if the pre- 
diction looks approximate. The not perfect agreement 
between theory and numerical calculation can be traced 
back to two major approximations. The first is the finite- 
ness of h which makes Eq. dl3) slightly inaccurate, while 
the second and the most important would be the ap- 
proximation related to consider the unperturbed states 
in Eq. (^. This last is in fact a double approximation, 
because it disregards the effect of the perturbation, but 
this is not so important as shown in Ref. |9[|, and the 
effect of the unperturbed decay which actually destroys 
the discrete levels picture we used. Nonetheless we think 
that our results are quite clear and to better show the re- 
lation between the peaks in the decay and the presence of 
nonlinear resonances disturbing the initial state we shall 
now resort to a graphical picture. 



V. A PHASE-SPACE REPRESENTATION: 
QUANTUM-CLASSICAL COMPARISON 

As we said before, in order to connect the quantum 
dynamics to the classical phase-space characteristics we 
must extend the concept of phase space to the quantum 
case. This can be done by using a phase-space represen- 
tation of the quantum state and among all the different 
possibilities we chose to use the Husimi representation, 
defined as: 



P{<1,P) 



dxaq^p{x)'ip{x)\\'^ , 



(17) 



where \aq^p) is a minimum indetermination state (coher- 
ent state) centered in {q,p). Using Eq. ( p^ we can ob- 
tain a phase representation of a quantum state in terms 
of the positive definite distribution p{q,p). In particu- 
lar we can calculate the Husimi distribution of the initial 



state which we chose to be an eigenstate of the unper- 
turbed well. This choice will allow us to obtain the cor- 
respondence we are looking for in a easy way, in fact in 
our approximation, namely, for small decay probability, 
we can safely disregard the dynamics inside the well for 
the times we explored. This means that the phase-space 
representation of the quantum state practically does not 
change during the time interval considered in Fig. |5|, and 
that in the comparison between quantum and classical 
phase space we can limit ourselves to deal with the ini- 
tial quantum distribution. 

In Fig. we show portraits of the classical phase space 
for increasing values of f2. For clarity we draw only the 
stochastic web and the nonlinear resonances island struc- 
tures, all the rest of the phase space being filled with reg- 
ular tori. As explained in Section pi due to the choice of 
the particular form of the time dependent perturbation, 
to the first order in e we have only resonances of the form 
i-Uq{E) — fl/n where ujq{E) is the frequency of the motion 
inside the unperturbed well. The frequency ujq{E) is a de- 
creasing function of the energy E for < -E < 2Vo, being 
equal to •\/4Vo for E = Q and vanishing for E = 2Vo which 
corresponds to the separatrix motion, see Eqs. (ra). This 
means that as we approach the separatrix we find reso- 
nances of larger order n. We can also easily realize that 
as the value of fJ is increased the nonlinear resonances 
move inside the phase space getting closer to the center 
of the well, and eventually disappearing when the rela- 
tion (ph cannot be fulfilled any more. In the meanwhile 
new resonances appear, moving out from the stochastic 
layer which is the region of the overlapping of all the 
infinite resonances of higher order n. 

In this motion towards the center of the well, the vari- 
ous resonances cross the region of the phase space which 
is occupied by the Husimi distribution of the initial quan- 
tum state (shaded area in Fig. 0b) and, as discussed in 
Section |^, we expect this to be related to the peaks of 
Fig. 1^. The region of parameters explored in Figs. ^ and 
is the same and from a first inspection we actually real- 
ize that the number of resonances crossing the shaded 
area corresponds to the number of peaks in the decay 
rate. This first result is already convincing but we can 
be more precise by singling out the phase-space portraits 
which correspond to the decay peaks. 
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FIG. 7. Classical phase-space portraits for differ- 
ent values of the driving frequency f2. The val- 
ues of the parameter are Vb = 0.048, e = 0.005, 
ft = 0.025. The value of f2, increasing from a) to j), is 
Q. = 1.4, 1.6, 1.7, 1.8, 1.9, 2.1, 2.2, 2.4, 2.5, 2.7, The 
shaded area in Fig. b) indicates the area under the Husimi 
distribution of the initial state, which in this case is the fourth 
eigenstate of the unperturbed well. 



This is done in Fig. |g where we show the classical 
phase-space structures and we indicate by the two con- 
tinuos lines the borders of the Husimi distribution of the 
initial quantum state. 
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The results of Section IV have been thus confirmed 
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FIG. 8. Classical phase-space portraits for different values 
of the driving frequency f2. The values of the parameter are 
Vo = 0.048, £ = 0.005, h = 0.025. The values of fi are 
chosen as the values corresponding to the peaks in Fig. Bp,, 
i.e., f7 = 1.73(a), 2.015(b), 2.25(c), 2.44(d). The two continu- 
ous circular curves indicate the outer and inner border of the 
Husimi distribution of the initial quantum state (the fourth 
eigenstate). The border is defined arbitrarily as the contour 
level of the distribution at 0.2 times the maximum height. 

The expected results is confirmed, the peaks in the 
quantum decay rate are related to the modification of the 
classical phase space under the quantum initial distribu- 
tion by a nonlinear resonance. The perturbation seems to 
be effective when a nonlinear resonance enter the exter- 
nal border of the initial state distribution and not when 
the nonlinear resonance passes over the center of the dis- 
tribution where the quantization torus is located. This 
is qualitatively in agreement with Eq. (nsh which pre- 
dicts that the perturbation appears when the nonlinear 
resonance is close to the semiclassical quantization torus, 
how much close being dictated by the quantum correc- 
tion which modifies Eq. (nsh with respect to Eq. (|5[). In 
fact the quantum correction to the forcing frequency fi is 
positive (the term "^j^ ' is negative for energy smaller 
than the barrier heights and thus the renormalized fre- 
quency is always larger than O) and this means that the 
decay peaks should appear for values of O smaller than 
the values for which the nonlinear resonances superpose 
to the quantization torus. This is precisely what happen: 
the peaks correspond to the approaching of the nonlinear 
resonances from the outside. 

This result comes along with a simple intuitive expla- 
nation: when a nonlinear resonance enters the external 
border of the initial state part of the initial distribution 
is moved outward by the islands structure, this produce 
a higher probability of tunneling across the barrier and 
thus an increase of the decay. This effect becomes less 
important once the nonlinear resonance penetrates inside 
the shaded region; until it disappears when the resonance 
is completely embedded in the central part of the distri- 



by the comparison of this section, making the picture of 
the role of the classical nonlinear resonances clearer and 
clearer. 



VI. CHAOS ASSISTED DECAY VERSUS CHAOS 
ASSISTED TUNNELING 

When we reviewed the numerical results we noticed 
that we found only enhancement of the decay compared 
to the unperturbed condition. This is not the case in 
the CAT, where the avoided crossings can produce both 
enhancement and decrease of tunneling, which can also 
vanish for particular values of the parameter |3[|. For 
this reason we wrote in the Introduction that the term 
assisted used in CAT is not really appropriate, but it 
seems that it could be better used in the present contest. 
To explain this different behavior we conjecture that this 
can be seen as a consequence of the continuity of the 
spectrum in our system. In fact the quenching of the 
tunneling is produced by the accidental degeneracy of 
the levels of the tunneling doublet, as seen in Fig. |^. 
This degeneracy is due to the modification induced by the 
crossing with the third level. In our case we do not have 
discrete levels, but a continuous density of states and thus 
the former picture simply does not apply. In other word 
for every value of the parameter the modification of the 
spectrum due to the presence of avoided crossings cannot 
lead to a complete degeneracy of the states and thus to a 
complete quenching of the decay. A similar situation and 
a graphical representation of this process can be found in 
a recent paper Q. On the other hand this continuous 
spectrum characteristic cannot completely rule out the 
possibility that chaos could produce also a decrease of 
the unperturbed decay rate, which we think could be 
present in some region of the parameters. 



VII. CONCLUSIONS 

A numerical calculation of the decay probability due to 
tunneling from a potential well showed that in presence 
of classical chaos the decay can be strongly enhanced 
and that this enhancement depends on the system pa- 
rameters in a resonant-like way. A qualitative inspection 
of the classical phase-space structure revealed a connec- 
tion between the peaks in the decay probability and the 
presence of classical nonlinear resonances in the region of 
the phase space occupied by the Husimi distribution of 
the initial state. This correspondence has been quanti- 
tatively explained using a semiclassical result which has 
been shown to be valid in the case of Chaos Assisted 
Tunneling. We can thus conclude that the enhancement 
of decay from a driven well can be explained by means of 
a perturbation of the region of the phase space occupied 
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by the initial state, perturbation produced by the emer- 
gence of nonhnear resonances. This is a direct connection 
between the modification of a purely quantum effect, the 
tunneling, and the classical phenomenon of destruction 
of the integrable dynamics which is at the basis of the 
chaotic behavior. 
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